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Nano-structuring impurity-doped crystals affects the phonon density of states and thereby mod¬ 
ifies the atomic dynamics induced by interaction with phonons. We propose the use of nano¬ 
structured materials in the form of powders or phononic bandgap crystals to enable or improve 
persistent spectral hole-burning and coherence for inhomogeneously broadened absorption lines in 
rare-earth-ion-doped crystals. This is crucial for applications such as ultra-precise radio-frequency 
spectrum analyzers and optical quantum memories. As an example, we discuss how phonon en¬ 
gineering can enable spectral hole burning in erbium-doped materials operating in the convenient 
telecommunication band, and present simulations for density of states of nano-sized powders and 
phononic crystals for the case of Y2Si05, a widely-used material in current quantum memory re¬ 
search. 


Rare-earth-ion-doped crystals have been studied for 
decades because of their unique spectroscopic properties 
that arise since their 4/-electrons do not partic¬ 
ipate in chemical bonding. At cryogenic temperatures 
they can offer narrow linewidths Hi] together with the 
possibility to spectrally tailor their broad inhomogeneous 
absorption lines [2]. These properties have led to many 
applications, including optical quantum memories EHa, 
signal processing [iniiiii, laser stabilization [HHH] , as 
well as ultra-precise radio frequency spectrum analyzers 
[isi [ig. Quantum memory implementations in rare- 
earth-ion-doped (REI-doped) crystals, such as those 
based on electromagnetically-induced transparency El), 
atomic frequency comb m, or controlled reversible 
inhomogeneous broadening [19], crucially rely on long 
coherence and spin-state lifetimes to achieve high 
efficiency and long storage time. 

Operating REI-doped crystals at low temperatures 
generally improves material properties. Yet, even at 
temperatures below 2 K, spin-lattice relaxation, i.e. 
thermalization of spins via interaction with phonons, 
still restricts lifetimes of spin states, reducing the ability 
to spectrally tailor the material. Eurthermore, by con¬ 
tributing to spectral diffusion [20] [21] and two-phonon 
elastic scattering processes [22], lattice vibrations limit 
coherence times. We propose the use of nano-structured 
materials to overcome these limitations. This is achieved 
by tailoring the phonon density of states [23] to restrict 
phonon processes, an approach not limited to REI-doped 
materials but also applicable to other impurity-doped 
materials such as color centers in diamond [24]. The 
result is an improved performance for all applications 
based on spectral hole burning, or that require long 
coherence times, by providing long-lived spin states while 
simultaneously suppressing phonon-driven decoherence. 


We note that the modification of population dynamics 
between REI crystal field levels in small powders, 
possibly related to phonon restriction, has been observed 
in the context of understanding luminescence dynamics 
[25l|26] - but not yet to improve spectral hole burning 
or optical coherence properties. 

In this paper, we propose a general procedure to 
identify the frequency range of detrimental lattice 
vibrations and restrict the relevant phonon processes in 
impurity-doped solids. As an example, we apply this 
approach to erbium-doped materials. Eirst, we recall the 
mechanisms causing detrimental electronic spin flips and 
discuss how the suppression of the direct phonon process 
influences spin-state lifetimes and spectral diffusion. 
Then, two methods for nano-structuring REI-doped 
crystals to suppress the direct phonon process affecting 
electronic spin states in Er^+:Y 2 Si 05 are analyzed: 
powder materials, which introduce a cutoff frequency 
in the phonon density of states, and phononic band 
gap crystals with a structure optimized for the targeted 
material and application. Numerical simulations for 
Er3+:Y2Si05 are presented. 

Relaxation of optically pumped spin populations back 
to thermal equilibrium by spin flips is caused by phonons 
primarily via three different processes [23 EH]: the di¬ 
rect process where phonons resonant with the spin tran¬ 
sition induce a spin flip, the Orbach process where two 
resonant phonons produce a spin flip via excitation to, 
and relaxation from, excited crystal field levels, and the 
off-resonant Raman process where two phonons lead to 
a spin flip via inelastic scattering. Levels involved in 
spectral hole burning, typically hyperfine or Zeeman lev¬ 
els, are mostly separated by MHz to GHz. The rates of 
phonon processes are proportional to the phonon density 
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FIG. 1. Magnetic field dependence of the spin-flip rate. As 
described by Equation the total rate R (solid red line) is 
composed of a constant term Ro (dotted blue line), a spin 
flip-flop term (green dashed line) and a direct phonon term 
(orange dash-dotted line), dominant above 1 T. We assume 
T = 3 K, ^ = 14, Ro = 0.1 kHz, ar> = 5 x 10~^ kHz/T^ and 
Qff = 2 kHz. 


of states in the material as discussed in m- Interest¬ 
ingly, the Debye approximation for the density of states 

p{uj) (X Vuj‘^ (1) 

with V the volume of the particle and uj the frequency, 
breaks down for small structures. In this regime, the den¬ 
sity of states becomes discrete [29] and acquires a cutoff 
frequency below which no vibrational modes are possible 
and thus no phonons can exist [23]. We predict that an 
ensemble of sufficiently small REI-doped crystals, or an 
appropriate REI phononic band gap material, will allow 
suppression of the direct phonon process responsible for 
relaxation between the spin states. 

As an example, we consider extending lifetimes of the 
ground-state Zeeman levels for persistent spectral hole 
burning and reducing spectral diffusion in erbium-doped 
crystals, particularly appealing for quantum memory 
applications due to the ^/i 5/2 A 3/2 transition of 

Er^+ ions around 1530 nm and long coherence times. 
Three important cases are Er^+:Y 2 Si 05 , for which 
optical coherence lifetimes of up to 4.4 ms have been 
observed [30], Er^+iLiNbOs, exhibiting a large oscillator 
strength accompanied by a very broad inhomogenous 
linewidth [31], and Er^+:KTi 0 P 04 , which has shown 
narrow homogeneous linewidths down to 1.6 kHz [32]. 
All Er-doped crystals so far eluded efficient quantum- 
state storage ED [ 33 H 35 ] since no efficient persistent 
hole burning involving electronic Zeeman levels [36] 
has been achieved. The only material that has shown 
long spin-state lifetimes is an Er-doped glass fiber 
ED, in which disorder broadens the spin transition and 
effectively weakens the coupling between Er^+ ions. This 
reduces the relaxation rate at low magnetic field. Intro¬ 
ducing controlled disorder in crystals is potentially an 


alternative approach to improve spin lifetimes however, 
disorder can be extremely detrimental for the coherence 
properties, as it is the case in impurity-doped glasses. 
Indeed, the latter are restricted in their applications 
since coherence times are intrinsically limited to a few ps 
or less by low frequency disorder mode dynamics that are 
present in any amorphous material [38] . The goal of our 
study is to propose a method to simultaneously achieve 
long coherence and spin lifetimes, which excludes, e.g., 
Er-doped fibers - at least given current knowledge. 


Phonons can limit the lifetime of persistent spectral 
holes, but other processes can also contribute, for exam¬ 
ple paramagnetic interactions of REI with neighbouring 
ions or with other impurities in the host. For erbium- 
doped materials such as Er^+:Y 2 Si 05 , Er^+iLiNbOs or 
Er^+:KTi 0 P 04 , the total spin-flip rate at temperatures 
below 4 K can be approximated by [20l[3ll|32] 


R{B, T) = Ro + sech^ 


V 2kT ) 


+ q^d g^B^ coth 


/ giiBB \ 
V 2kT ) 


( 2 ) 


The second term represents the average rate of mutual 
spin flip-flops, i.e. exchange of spin states between two 
Er ions, with a constant, g the factor of the ions, 
Pb the Bohr magneton and k the Boltzmann constant. 
The last term corresponds to the direct phonon process 
and is proportional to for gpBB < 2/cT, where the 
parameter q^d characterizes the strength of the phonon 
coupling ED- The constant Rq includes all residual 
relaxation processes. The Raman and Orbach rates are 
not considered here since they are strongly suppressed at 
low temperature [28], particularly when the crystal-field 
splitting is much larger than kT as is the case for 
the materials considered here. The magnetic field 
dependence of the spin-lattice relaxation is illustrated 
in Figure [^ which shows the total rate and its different 
components. Values of the parameters vary for the three 
materials cited above. The Rq term is negligible in 
EriLiNbOs [SI], whereas the spin flip-flop rate is very 
small for Er:Y 2 Si 05 [20]. Er:KTi 0 P 04 has unusually 
small Rq and weak spin flip-flop process, so that only 
the direct phonon process significantly contributes to 
the field-dependent spectral diffusion, even at low field 
strengths [32]. However, in all cases the direct phonon 
process becomes dominant for magnetic field strengths 
of 1 T and higher. 


From the above discussion, the optimum operating 
point would seem to be at magnetic fields below 0.5 T; 
however, to reduce spectral diffusion and improve coher¬ 
ence properties, erbium-ion-doped-crystals often require 
magnetic fields of a few Tesla. It has been shown that 
the maximum linewidth due to spectral diffusion can be 
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FIG. 2. Effect of the direct phonon process suppression on the 
effective homogeneous linewidth. As described by Equation 
Feff after t = 10 /xs is shown as a function of the magnetic 
field strength in the presence (blue dashed line) and in the 
absence (red solid line) of the direct phonon process. We 
assume the same parameters as in Eig. and, in addition, 
r^ax = 5 X 10 ^° Hz, Fo = 10 kHz. 


described by m 


rsD(B,r) = r„,,sech" ’ (3) 

with Fmax the full-width at half-maximum frequency 
broadening caused by magnetic dipole-dipole interactions 
between Er ions. One can see that Fsd decreases with 
magnetic field and spectral diffusion can thus be sup¬ 
pressed at high fields due to the increased magnetic or¬ 
der in the system. This leads to narrow homogeneous 
linewidths at high fields. However, high fields lead to in¬ 
creased spin-lattice relaxation (see Figure [^. In addition 
to being detrimental to persistent spectral hole burning, 
this results in faster broadening of the spectral line to its 
maximum value, as shown in the effective linewidth after 
a time delay t m 

Tesit) = To + Tsd (1 - e-«‘) , (4) 

with Fq the homogeneous linewidth without spectral 
diffusion. Therefore, removing the direct phonon process 
will not only improve persistent spectral hole burning, 
but also inhibit spectral diffusion, especially at magnetic 
fields around 1 T (see Figure]^. 

We now identify the range of frequencies at which 
phonons must be suppressed. In our example of erbium- 
doped materials, the goal is to eliminate phonons below 
the frequency corresponding to the maximum Zeeman 
splitting under an external magnetic field of around 
1 T. In Er:Y 2 Si 05 the ^-factor along the 6-axis is 
13.6 [39], leading to a splitting of 190 GHz/T. As we 
will describe next, inhibiting phonons at such high 
frequency requires small structures. For other types of 
paramagnetic centers or nuclear hyperfine states, the 


energy splitting of the spin states is usually smaller, 
so that larger structures can inhibit the detrimental 
phonon processes. Furthermore, it is possible to suppress 
phonons resonant with transitions between crystal-field 
levels, depending on the energy level splitting. In 
consequence, two-phonon processes such as Raman or 
Orbach processes can be suppressed. This would also 
lead to improved coherence lifetimes for materials in 
which two-phonon elastic scattering is responsible for 
dephasing [22] . More generally, suppressing phonons will 
enable to increase the range of working temperatures 
for all impurity-doped materials that have properties 
limited by phonon processes. 

We first discuss the case of small crystals. To sup¬ 
press phonon processes we use the fact that the density 
of states of a particle features a cutoff at 

c 

^min Vmin j i 
ird 

with c the sound velocity, d the diameter of the 
crystal, and r]min a numerical constant that is 2.05 
for a spherical particle [40] . Below this cutoff no 
phonons are allowed. In the region above the cutoff, the 
phonon density of states exhibits gaps before it transi¬ 
tions into a continuous function described by Equation [^ 

To simulate the phonon density of states for a spherical 
particle, the eigenmodes must first be calculated. They 
can be separated into spheroidal-transverse, spheroidal- 
longitudinal, and torsional modes. The characteristic 
equations needed to calculate their frequencies have been 
described in a quantized form by Takagahara m- Each 
eigenmode is identified using three quantum numbers: 
the radial quantum number j, the angular momentum 
/, and cr, which specifies the type of mode. Using those 
eigenmodes, the phonon density of states can be written 
as [42] 


P(^) — ^(2^ + 1) 




Auj 




(6) 


Each eigenmode is 2/ + 1 times degenerate. We assume 
that the density of states is composed of Lorentzian 
shaped lines at the frequency of the respective eigenmode 
cjU and linewidth Auj. In the work of Meltzer and Hong 
[12 , the Lorentzian width was assumed to increase with 
the frequency squared, an empirical assumption that 
described their particular experimental data adequately. 
However, since there is no theoretical framework to 
support that assumption, the width is assumed to be 
constant here. 


In the case of Er^+:Y 2 Si 05 around 1 T, we estimated 
from Eq. [^ that nanoparticles of 12 nm diameter are 
suitable to suppress the direct phonon process. The 
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FIG. 3. Simulated phonon density of states in a Y 2 Si 05 
nanoparticle with a diameter of 12 nm. Each eigenmode of 
the spherical particle is represented by a Lorentzian with a 
width of 1 GHz. 


phonon density of states for a Y 2 Si 05 nanoparticle of 
12 nm diameter, calculated according to Equation 
is shown in Figure Phonon processes at frequencies 
below 200 GHz are fully restricted. This assumes that 
there is no phonon propagation between the particles. 
Simulating the effect of touching particles on phonon 
propagation is not trivial and subject further theoretical 
and experimental research. However, for a lightly packed 
powder immersed in cold helium gas we expect proper¬ 
ties approaching those of single, isolated particles. Note 
that particles separated by electrostatic interaction can 
be obtained by coating the powder with the appropriate 
surfactant or by specific syntheses m- 

In powders, averaging of anisotropic properties over 
the random distribution of crystallite orientations 
generally causes issues in applications such as quantum 
memories and spectroscopy. However, if the material 
has a dielectric or magnetic anisotropy, the crystals can 
be preferentially oriented using external fields [44l |45] . 
This allows the control of parameters, such as Zeeman 
splittings and the direction of the ion’s transition 
dipole moment. Scattering of laser light by powder 
samples can be reduced or even eliminated if the size 
of the crystals is smaller than the wavelength of the 
interacting light, since the powder can be treated as a 
bulk material with an effective (average) refractive index. 

Crystal powders can be produced using several 
approaches and spectroscopic properties comparable 
to bulk crystals have been observed in some cases 
[46] . A thorough review of different growth methods 
can be found in [47] and references therein. Growing 
nano/microcrystals has the advantage of being relatively 
inexpensive. An alternative way to produce nanocrystals 
is to grind a bulk crystal into powder using techniques 
such as ball milling. With this approach, sizes on 
the order of 10 nm are achievable by using very small 
grinding balls IlHl ES]. However, reducing dimensions 
below around 10 nm may result in additional (detrimen¬ 



FIG. 4. (a) Unit cell of the proposed phononic crystal struc¬ 
ture. The radius of the island is 9.5 nm, the width of the arms 
is 5 nm, and their length 5 nm. This structure is 7 nm thick, 
(b) Band structure of the phononic crystal in (a). The band 
gap (shaded grey) around the center frequency of 178 GHz is 
clearly visible, (c) Reciprocal lattice [59] showing the points 
traversed in the band structure. 

tal) decoherence [42] [50] and lifetime limiting effects [51] . 

A second possibility to avoid detrimental phonons is to 
create a phononic band gap material [52H54] engineered 
to produce a gap in the density of states around the 
undesired frequencies. Such structures are employed 
in the field of optomechanics [55] |56], and are also 
used for thermal phonon engineering i.e. for thermal 
metamaterials and thermocrystals EH- By changing the 
geometry and overall size of phononic crystal structures, 
it is possible to optimize the band gaps for a specific 
application [58] . 

Phononic crystal materials have the advantage of 
employing an oriented single crystal. In addition, the 
problem of contact between nanoparticles would be 
avoided. The fabrication of phononic band gap materials 
has already been explored for materials such as silicon 
or gallium arsenide m and phononic structures with 
a bandgap at 33 GHz and dimensions around 100 
nm have been demonstrated m- To create phononic 
crystal structures in currently-used host materials such 
as Y2Si05, Y3AI5O12 or LiNbOs, Ihin films of those 
materials would be a prerequisite. Thin films of LiNbOs 
are commercially available and there have been demon¬ 
strations of thin single-crystal film fabrication for other 
REI host materials [62] [63]. Finally, structuring these 
materials with nanometer precision has been shown to 
be feasible using focused ion beam milling [64] [65] . 

To illustrate this approach for our example of EriYSiO, 
we designed and optimized a phononic structure using 
finite element simulations (COMSOL). The structure 
is inspired by designs of Safavi-Naeini et al.[66]. We 
engineered our structure such that it produces a band 
gap around 178 GHz with a width of 3.5 GHz. Figure 
[^ shows the proposed geometry of the phononic crystal 
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unit cell together with its band structure. Even though 
it is still challenging, the rapid progress in the field of 
nano-fabrication together with the discovery of helium 
focused ion beam milling allows structures with these 
dimensions to be fabricated out of a thin film [671 EHl- 
Through additional etching it is possible to remove 
part of the substrate below the phononic crystal such 
that the structure suspended and only connected to the 
substrate at its borders. This ensures complete phonon 
suppression. 

In summary, we propose to restrict phonon processes 
in REI-doped materials through nano-tailoring in order 
to improve persistent spectral hole burning and reduce 
spectral diffusion. We briefly discussed two methods to¬ 
wards this end - the use of nano-crystals and engineered 
phononic band gap structures - and described how our 
concept can be applied to well-studied Er-doped crys¬ 
tals. In turn, this would allow realizing optical quan¬ 
tum memories operating at convenient telecommunica¬ 
tion wavelengths and increase the performance of other 
applications that are based on persistent spectral hole 
burning or that rely on narrow homogeneous line widths. 
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